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Background: Following the 2007 precise measurements of monopole strengths in tin isotopes, there has been a continuous 
theoretical effort to obtain a precise description of the experimental results. Up to now, there is no satisfactory explanation 
of why the tin nuclei appear to be significantly softer than ^"^Pb. 

Purpose: We determine the influence of finite-range and separable pairing interactions on monopole strength functions in 
semi-magic nuclei. 

Methods: We employ self-consistently the Quasiparticle Random Phase Approximation on top of spherical Hartree-Fock- 
Bogolyubov solutions. We use the Arnoldi method to solve the linear-response problem with pairing. 

Results: We found that the difference between centroids of Giant Monopole Resonances measured in lead and tin (about 
1 MeV) always turns out to be overestimated by about 100%. We also found that the volume incompressibility, obtained by 
adjusting the liquid-drop expression to microscopic results, is significantly larger than the infinite-matter incompressibility. 

Conclusions: The zero-range and separable pairing forces cannot induce modifications of monopole strength functions in tin 
to match experimental data. 

PACS numbers: 21.60.Jz, 21.10.Pc, 21.10.Fe 



I. INTRODUCTION 

The incompressibility of infinite nuclear matter as well 
as of finite nuclei has been studied in a number of the- 
oretical papers and reviews. In the classic review by 
Blaizot |l| the connection between the finite-nucleus in- 
compressibility and centroid of the Giant Monopole Res- 
onance (GMR) was shown. This relation allows us to 
study incompressibility of nuclei through microscopic cal- 
culations of the monopole excitation spectra. It also 
brings us the possibility to directly compare theoretical 
results with experimental data. For examples, see the 
measurements presented in Rcfs. 043- 

In Ref. 0, it was shown that the self-consistent mod- 
els that succeed in reproducing the GMR energy in the 
doubly-magic nucleus ^"^^Pb systematically overestimate 
the GMR energies in the tin isotopes. In spite of many 
studies related to the isospin , surface , and pair- 
ing [lol - [l6j influence on the nuclear incompressibility, to 
date there is no theoretical explanation of the question 
"Why is tin so soft?" [1, For an excellent recent re- 
view of the subject matter we refer the reader to Ref. 3. 

Studies in Rcfs. [l^, [l^ were restricted to the ef- 
fect of zero-range pairing interaction. In the present 
paper we focus on a different kind of pairing force, 
namely, we implement the finite-range, fully separable, 
translationall y in variant pairing interaction of the Gaus- 
sian form |18l - [2(J | , together with the general phenomeno- 
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logical quasilocal energy density functional in the ph- 
channel [2l|. We have performed calculations for all 
particle-bound semi-magic nuclei starting from Z = % 
or TV = 8, up to Z = 82 or = 126. The ground- 
state properties were explored within the Hartree-Fock- 
Bogolyubov (HFB) method, whereas the monopole exci- 
tations were calculated by using the Quasiparticle Ran- 
dom Phase Approximation (QRPA) within the Arnoldi 
iteration scheme [22j. For the numerical solutions, we 
used an extended version of the code HOSPHE [23j . 

The paper is organized as follows. In Sees. Ulland IIIII 
we briefly outline the Arnoldi method to solve the QRPA 
equations and present the separable pairing interaction, 
respectively. In Sec. IIVI we discuss the nuclear incom- 
pressibility, including its theoretical description, deflni- 
tions in finite and infinite nuclear matter, and relations 
to monopole resonances. Then, our results are shown and 
discussed in Sec.lVland conclusions are given in Sec. IVII 



II. QRPA METHOD 

In the present study, we solve the QRPA equations 
by u sing the iterative Arnoldi method, implemented in 
Ref. [221 ■ It provides us with an extremely efficient and 
fast way to solve the QRPA equations. The QRPA equa- 
tions are well known [23, and have been recently re- 
viewed in the context of the finite amplitude method [26} . 
Therefore, here we only give a brief resume of basic equa- 
tions, by presenting their particularly useful and compact 
form. 

Basic dynamical variables of the QRPA method are 
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given by the generalized density matrix TZ, 

corresponding to mean- field Hamiltonian H = d£/dTZ, 
'h-X A 



n = 



A+ -h* + x r 



(2) 



The standard HFB equations that define amplitudes U 
and V read 



h-\ 
A+ 



A 



U V* 
V u* 



U V* 

V u* 



E 
-E 



. (3) 

where the diagonal matrix E contains positive quasipar- 
ticle energies. Then the quasiparticle (%) and quasihole 
{lp) states are given by columns of eigen- vectors: 



X ■= 



(4) 



that is, 

Hip = -ipE, Hx = XE. (5) 

The vibrational time-dependent HFB state 

|*(t)) = |vI/)+e-*|*), (6) 

where |^) is a small- amplitude correction, leads to the 
time-dependent density matrix, 

7^(^) = 7^ + c*"*-^ + e-''^*7t+ (7) 

and time-dependent mean field H(t), 

nit) ^ n + e"^'n + e-"^'n+ . (8) 

After a linearization of fields in the time-dependent 
Hamiltonian, one obtains the QRPA equations in a sim- 
ple form. 



hion = [H,n] + [H,n]. 



(9) 



In this approach, states in Eq. (jG]) play a role of Kohn- 
Sham-like wave functions, which serve the purpose of 
generating generalized density matrices TZ{t) only. Nei- 
ther represents a correct ground state of the system 
nor \^) represents that of an excited vibrational state. 
However, the amplitude TZ, which constitutes the fun- 
damental degree of freedom of the QRPA method, docs 
represent a fair approximation to the transition density 
matrix between both states of the system. It then allows 
for calculating matrix elements of arbitrary one-body op- 
erators between the ground state and vibrational state, 
which is the primary goal of the QRPA approach. 

Equation ^ constitutes the base for our solution of 
the QRPA equations in terms of the iterative Arnoldi 



method. Indeed, since the mean-field amplitude H de- 
pends linearly on the density amplitude TZ, Eq. ^ con- 
stitutes an eigen-equation determining TZ and hui. How- 
ever, the matrix to be diagonalized, that is the QRPA 
matrix, does not have to be explicitly determined. To 
obtain the entire QRPA strength function, it is enough 
to start from a pivot amplitude and repeatedly act on it 
with the expression on the right-hand side |22| . In each 
iteration, one only has to calculate the mean-field ampli- 
tude Ti corresponding to the current density amplitude 
TZ, which is an easy task. The pivot can be freely chosen 
to optimally suit the calculation. It can for example be 
random, a QRPA eigen-phonon or be constructed from 
an external field. In this work we construct the pivot 
from the monopole transition operator. This approach is 
fundamentally different than that used within the FAM 
of Ref. [2^, where an external field is used throughout 
the calculation and Eq. ^ has to be iterated for all val- 
ues of frequencies uj. 

Since both stationary (7?.^ = TZ) and time-dependent, 
{TZ^{t) = TZ{t)) density matrices are projective, the 
QRPA amplitude TZ has vanishing matrix elements be- 
tween the quasihole and between the quasiparticle states, 
that is, 



ip+TZip = x+^X = 0. 



(10) 



Therefore, TZ is solely defined through the antisymmetric 
amplitude matrices Z and Z'^ defined as 

z = -z^ = x^np, 
z'+ = -z'* = p+izx- (11) 

Explicitly, amplitudes Z and Z''^ read 

Z = U+pV* + U+kU* + V+k'+V* -V+p^U*, 



Z'+ = V' pU + V' kV + k'+U - U^ p'V. (12) 

Within such a formalism, the QRPA equations © can 
be expressed as 



-hwZ = EZ + ZE + W, 
hujZ'+ = EZ'+ + Z'+E + W'+ , 



(13) 



where the field amplitudes W and W'^ are defined as 

W'+ = -W'* = ip+iix, (14) 

or explicitly, 

W U+hV* + U+AU* + V+A'+V* -V+h^U*, 
W'+ = V'^hU + V'^AV + U'^A'+U -U'^h^V. (15) 

We can also invert Eq. ([T^ and obtain transition den- 
sities p, k, and k'~^ expressed in terms of amplitudes Z 
and that is, 



p = UZV^ + v*z'+u+. 



k = UZU' + V*Z'+V+ 
k'+ = VZV^ + u*z'+u^ 



(16) 
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Finally, we can reduce the above QRPA formalism to 
spherical symmetry used in the present study. Then, 
the vibrating amplitude of Eq. ([6|) has good angular- 
momentum quantum numbers JM, that is, |^') = [vj/'^*^) 
and hence all the QRPA amplitudes pertain to the given 
preselected channel JM, while the ground state \^) is 
spherical. As a consequence, as dictated by the angular- 
momentum algebra, only specific spherical single-particle 
states are coupled by the QRPA amplitudes, which can 
be expressed through the Wigncr-Eckart theorem and re- 
duced matrix elements as 

1 „„„ . , , , ^^^^ 



X 



JM 

ajm.a' j' m' 



V27+T 



^j'Tn'JA/(V'ajl 



where X stands for amplitudes p or h, and 

^ajm,a' j'rn' ^ ^ 2 J + jmj' m' {''PajWX Wlpa'j'), (18) 

>' + J\ 



^ ajm,a'j'm' " ^ 2 J + f^^3"'i'™' ^^"^^ 1 1 V^a'i') (19) 



where X stands for amplitudes k, A, Z, or W. In these 
expressions, we have used the standard quantum num- 
bers ajm of spherical single-particle states. 

Spurious QRPA mode appears in the 0+ QRPA cal- 
culations. In a self-consistent full QRPA diagonaliza- 
tion, the spurious mode decouples from the physical 
QRPA modes and appears at zero energy. In the Arnoldi 
method, this separation does not happen unless we make 
the full Arnoldi diagonalization, which usually is not fea- 
sible. 

To prevent the mixing of physical QRPA excitations 
with the spurious 0^ mode, before the Arnoldi iteration 
we create the spurious-mode QRPA amplitudes and its 
associated conjugate-state (boost-mode) QRPA ampli- 
tudes. The spurious 0^ mode amplitudes follow from 
the particle number operator and have the form. 



pOO ^ U+V*, P' 



J+OO 



= V U. 



(20) 



The 0+ boost mode is generated by making an additional 
HFB calculation whose chemical potentials A,- and aver- 
age particle numbers are slightly shifted from the ground 
state values, producing a perturbed state IIIFB2). The 
boost-mode amplitudes are calculated by using Thouless 
theorem as, 



R 



w 



n+00 
ajm^a' j' m' 



(HFB2|a+,„a+^.,„,,|HFB) 
(HFB2IHFB) 

(HFB|aa'j'm'aQj,„|HFB2) 



(HFBIHFB2 



vu- 



(21) 



(22) 



where we used the standard transformation matrices 
from one quasiparticle basis to another p3 |. 



V = U^V2 + V^U2, 



u 



U+U2 + V+V2. 



(23) 
(24) 



Gram-Schmidt orthogonalization is used to keep dur- 
ing the Arnoldi iteration the Krylov-space basis vectors 
orthogonal to the spurious and boost modes, that is, each 
Krylov-space basis vector is orthogonalized against P and 
R. The orthogonalization procedure is described in detail 
in Rcf. [22|. For the semi-magic nuclei considered here, 
we only vary the particle number of the nucleon species 
that has non- vanishing pairing correlations. 



III. SEPARABLE PAIRING INTERACTION 

The separable finite-range pairing interaction for neu- 
trons {t ~ n) and protons (r = p) that we use in this 
study is defined as [l9| 

Vr{risi,r2S2;r[s[,r'2S2) 

- -Gr5{R - R')P{r)P{r')^{l ~ P„), (25) 

where R ~ (ri + r2)/2 denotes the centre of mass coor- 
dinate, r = ri — r2 is the relative coordinate, r = |r|. 
Pa is the standard spin-exchange operator, and function 
P{r) is a sum of m Gaussian terms. 



m ^ (47ra^)3/2 



(26) 



Coupling constants Gt define the pairing strengths for 
neutrons and protons. 

For such a pairing interaction, the pairing energy ac- 
quires a fully separable form, which in spherical symme- 
try reads 



pair 



1 



Gr 



NJt 



f+Jl 



and depends on the reduced matrix elements of the 
pairing densities k^- and k'^ between the single-particle 
wave functions ip^^r) for fi denoting the set of spheri- 
cal harmonic-oscillator quantum numbers Uf^lf^j^. The 



interaction matrix elements 1/^'^are defined as 



J" 



= y'(47r)(2J + l)(2j^ + l)(2j, + l)<j i i 



X M. 



NJnO 



2I/4 /7ri/2(2n + l)! 1 



1 



1 + aH^ 



63/2 y 2(2"n!)2 m ^ (47raf)3/2 

3/2/1 -a^fe^s 

^ ' ' (28) 



1 + a?62 



where 2n = 2n„ + L + + - 2N 



J ]\,rNXnO 



are the standard Talmi-Moshinski coefficients [27|, and 
b = ^Jrrujjjh, denotes the harmonic-oscillator constant. 
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IV. NUCLEAR INCOMPRESSIBILITY 

The isoscalar incompressibility of infinite nuclear mat- 
ter is defined by the well-known formula [l| 



9p' 



dp2 



AJ p=p„ 



(29) 



where p„,„ is the saturation density of nuclear matter. 
Of course, K^o cannot be directly measured; however, by 
using Eq. it can be calculated from theoretical equa- 
tion of state E{p) or it can be indirectly estimated from 
measurements of monopolc excitations of finite nuclei. 

The incompressibility of finite nucleus, K^, is defined 
by its scaling-model relation [2^ to the centroid of the 
giant monopole resonance (GMR), £^gmr, as 



GMR — 



m{r'^) 



(30) 



where (r^) is the average square radius of the nucleus. 
Eq. (f30| is derived under the assumption that most of the 
monopole strength is concentrated within one dominant 
peak, see Ref. The centroid of the GMR can be 

extracted from its strength function as the ratio of the 
first and zero moments, that is. 



Egmk ~ — -• 
mo 



(31) 



There exist several alternative ways to extract -Bqmr 
through differ ent mome nts of the strength function, such 
as i?GMR = \/mi/m_i or Eqmr = yjm^/mi. However, 
they are more sensitive to details of the strength function 
and thus less appropriate for studies of the incompress- 
ibility. 

In analogy to the Weizsacker formula for the nuclear 
masses, one can introduce [l| a similar relation for nuclear 
incompressibilities , 



^4/3 



(32) 



Similarly as in the liquid-drop model, we refer to Kvi 
Ks, Kt, Ks.t, and Kc as the volume, surface, sym- 
metry, surface-symmetry, and Coulomb incompressibility 
parameters, respectively. By adjusting these parameters 
to the incompressibilities Ka, calculated in finite nuclei 
from Eqs. ([50)1 and ([5T|) . we can obtain an estimate of 
the infinite-matter incompressibility as Koa — Ky. 



V. RESULTS 

In our study we performed a set of calculations for 
semi-magic nuclei starting from Z = 8 or N = 8 and end- 
ing with Z = 82 or iV = 126. The ground states proper- 
ties were calculated within the HFB method by using the 



code HOSPHE [23[, whereas the monopole strength func- 
tions were obtained by implementing in the same code 
the QRPA method within the Arnoldi iterative method. 

We decided to use two different Skyrme functionals ~ 
SLy4 [111 and UNEDFO Both of them were tuned 

(among other observables) to reproduce the main prop- 
erties of the infinite nuclear matter. In particular, they 
correspond to the same value of nuclear incompressibility 
(p9)) of Koo = 230 MeV and differ in their values of the 
effective mass of m*/m = 0.70 and 1.11 for SLy4 and 
UNEDFO, respectively. 

The present study is focused on comparing incom- 
pressibilities obtained with two different pairing interac- 
tions, namely, the standard zero-range force, Vr{r,r') 
— Vbr<5(r — r'), and separable force presented in Sec 
To make the comparison meaningful, we adjusted the 
strength parameters, Gr and Vr, so as to obtain for 
both forces very similar neutron (proton) pairing gaps 
in Z = 50 isotopes {N = 50 isotoncs). The resulting 
gaps roughly correspond to the experimental odd-even 
mass staggering along the Z = 50 and = 50 chains 
of nuclei. Theoretical pairing gaps. A, 
determined as in Ref. [3l|, namely. 



and Ap, were 



A, 



Trpr 



(33) 



where TtA = J2k^kk and TvA = Efc>o^fefe- Fo'" 
separable pairing, in Eq. (j26p we used only one Gaussian 
term with cii = 0.66 fm. 

In this way, in the calculations we used the separable- 
force strength parameters of G„ 631 and 473 MeV fm'^ 
(Gp = 647 and 521 MeVfm^) for the SLy4 and UNEDFO 
functionals, respectively, and similarly, for the zero-range 
force: U„ = 195 and 126 MeV fm^ {Vp = 221 and 
157MeVfm'^). All calculated neutron and proton pairing 
gaps are shown in Figs. [T] and [2l respectively. One can 
see that the results obtained for both pairing forces are 
fairly similar. The HFB iterations were carried out using 
a linear mixing of densities from the current and pre- 
vious iteration defined by a constant mixing parameter 
[2^. With this recipe, for some of the nuclei, the HFB 
iterations did not end in converged solutions. Such cases 
were excluded from the analysis of pairing properties and 
the subsequent QRPA calculations. 

We note here that no energy cut-off is needed for calcu- 
lations using the separable force, and thus in our calcula- 
tions the entire harmonic-oscillator basis up to Nq = 20 
shells was used. On the other hand, for the zero-range 
force we used the cut-off e nerg y of 60 MeV applied within 
the two-basis method [3^ . l33l | 

First; in Fig. [3] we illustrate the high reliability of the 
Arnoldi method in determining the key factors of our 
analysis, namely, the ratios of moments of the monopole 
strength functions. To obtain a perfectly stable result, 
only about 70 Arnoldi iterations suffice. In this way, the 
QRPA result is achieved within the CPU time that is of 
the same order as that needed to obtain a converged HFB 
ground state. Note that the Arnoldi iteration conserves 
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FIG. 1: (Color online) Neutron pairing gaps in the Z = 8, 
20, 28, 50, and 82 isotopes (see the legend shown in Fig. [5]). 
Upper and lower panels show results obtained for the SLy4 
and UNEDFO functionals, respectively. Left and right panels 
show results obtained for the zero-range and separable pair- 
ing, respectively. 
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FIG. 2: (Color online) Same as in Fig. [T] but for the proton 
gaps in the N = 8, 20, 28, 50, 82, and 126 isotones (see the 
legend shown in Fig. [B} . 



all odd moments, so during the iteration, the moment mi 
does not change; thus the convergence of mi /mo simply 
illustrates the convergence of mp alone. 

In Fig. 2] we compare our QRPA results with raw ex- 
perimental data obtained in Ref. In this work, a 
Lorentzian fit to data was performed in the region of 
energies of 10.5-20.5 MeV, and the experimental values 
of mi/niQ were determined from the corresponding fit- 
ted curve (its moments were calculated for energies from 
zero to infinity). In determining our theoretical values 
of mi /mo, we also perform the integration in the entire 
energy domain. We have checked that the integration of 



FIG. 3: Convergence of the ratio of first and zero moments 
mi/mo calculated in ^^^Sn as a function of the number of 
Arnoldi iterations. 



theoretical curves in the fixed region of 10.5-20.5 MeV 
does not bring meaningful results, because, in the wide 
region of masses studied here, the GMR peaks move too 
much, and extend beyond the above narrow range of en- 
ergies. Our QRPA strength functions were obtained from 
the discrete Arnoldi strength distributions by using the 
smoothing methods explained in Ref. (2^. We also note 
that in our QRPA calculations, the high-energy shoulder 
of the strength fimction is not obtained, cf. discussion in 
Ref. 0. 
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FIG. 4: (Color online) The QRPA monopole strength function 
in ^^^Sn (solid line) compared to raw experimental data 0] 
and Lorentzian fit to data (dashed line) performed in the re- 
gion of energies of 10.5-20.5 MeV 

Figs.[5]and|6]present the overview of all obtained finite- 
nucleus incompressibilities Kj^, Eqs. ([50)1 and (|3ip . calcu- 
lated along the isotopic and isotonic chains, respectively. 
One can see that for both Skyrme functionals, SLy4 and 
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UNEDFO, values corresponding to the zero-range (full 
symbols) and separable (open symbols) pairing forces are 
very similar. 



course, we can expect that using Skyrme parametriza- 
tions tuned to higher (lower) values of i^oo may lead to 
uniformly higher (lower) values of Ka- 
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FIG. 5: (Color online) Incompressibility Ka calculated for the 
isotopic chains of semimagic nuclei with Z — 8, 20, 28, 50, and 
82. Left and right panels show results obtained for the SLy4 
and UNEDFO functionals, respectively. Full (empty) symbols 
correspond to the zero-range (separable) pairing force. 



FIG. 7: (Color online) Incompressibility Ka calculated for 
chains of the Z — 50 and 82 isotopes. Results obtained by 
using the separable (squares) and zero-range (circles) pairing 
with the UNEDFO functional are compared to the available 
experimental data 
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FIG. 6: (Color online) Same as in Fig. [5l but for the isotonic 
chains with N = 8, 20, 28, 50, 82, and 126. 

To see effects of the pairing interaction in more de- 
tail, we focus on the results obtained for chains of tin 
and lead isotopes. In Figs. [7] and [5] we compare the- 
oretical results with the experimental data for ^°^Pb 
and ^^^~^^^Sn, taken from Refs. ^ comparison 

of the two types of pairing interactions, and two dif- 
ferent Skyrme functionals, leads to the conclusion that 
the calculated incompressibilities Ka depend on the in- 
teractions in the particle-particle channel as well as the 
particle-hole channel of the two Skyrme functionals used 
in our study - SLy4 and UNEDFO - only weakly. Of 
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FIG. 8: (Color online) Same as in Fig.0 but for the UNEDFO 
(squares) and SLy4 (circles) functionals and separable pairing 
force. 

To check a weak dependence of Ka on the inten- 
sity of pairing correlations, we have repeated the cal- 
culations by using values of neutron pairing strengths 
varied in a wide range, Gn = 631 ± ISOMeVfm'^ and 
Vn = lOSztSOMeVfrn-^. Such variations induce very large 
changes of neutron pairing gaps, shown in Fig.[9l the ones 
that are certainly beyond any reasonable range of un- 
certainties related to adjustments of pairing strengths to 
data. In Figs.fTUlandlTTl we show the influence of the var- 
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ied pairing strengths on the calculated incompressibilities 
Ka. We see clearly that even such large variations can- 
not induce changes compatible with discrepancies with 
experimental data. 
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FIG. 9: (Color online) Neutron pairing gaps calculated in tin 
isotopes for low (triangles), central (squares), and high (cir- 
cles) values of pairing strength parameters given in captions 
of Figs. [TOl and [TT] 
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FIG. 10: (Color online) Incompressibility Ka calculated 
for the SLy4 functional and separable pairing force in tin 
(squares) and lead (circles) isotopes compared to the available 
experimental data. Theoretical results are plotted together 
with uncertainties pertaining to variations of the neutron 
strength parameter in the range of G„ = 631 ± 150MeVfm'^ 



To illustrate the effect of isospin asymmetry, in Figs.fTOl 
and [11] we plotted the results as functions of N/Z, 
whereby ""^^^Sn and ^°^Pb arc located at almost the same 
point of the abscissa. These figures clearly show that 
the discrepancies with data are probably not related to 
the isospin dependence of Ka- Indeed, for both types of 
pairing, in the region of 1.0 < N/Z < 1.6, the results 
obtained for tin and lead isotopes roughly follow each 
other. 

Finally, to illustrate the fact that nuclear radii are 
fairly robust and cannot significantly influence the val- 
ues of Ka, determined from Eqs. ([50)1 and ([5T|) . we show 
values of toi/toq alone in Figs. [T^ and [T51 We see that for 
both types of pairing, in tin and lead the calculated values 
of mi /mo overestimate and underestimate the measured 
ones by 0.6-0.8 and 0.4 MeV, respectively. Exactly the 
same pattern was obtained within the relativistic nuclear 
energy density functionals studied in Ref . [l^ , where the 
corresponding discrepancies were equal to 0.8-1.0 and 
0.2 MeV. We also note that this comparison directly re- 
lates calculations to data, without using the intermediate 
and model-dependent definition of Ka- 

To conclude our analysis, we have also performed ad- 
justments of the liquid-drop formula p2p to our micro- 
scopically calculated values of Ka- The obtained param- 
eters are collected in Table ID We see that the liquid-drop 
formula is able to provide an excellent description of the 
QRPA results, with average deviations of the order of 
5 MeV, that is, about 3% of the typical value of Ka- 
Similarly the values of the volume incompressibility Ky 
are determined to about 2% of precision. The least pre- 
cisely determined liquid-drop parameter is the surface- 




N/Z 

FIG. 11: (Color online) Same as in Fig. 1101 but for the zero- 
range pairing force and uncertainties pertaining to variations 
of the neutron strength parameter in the range of Vn = 195 ± 
30 MeVfm^ 



symmetry incompressibility Ks^t, estimated up to 25% 
of precision. We also note that, within the fit precision, 
the volume parameter Ky averaged over both function- 
als and both pairing forces equals to 254±5 MeV, which 
is significantly higher than the corresponding infinite- 
matter incompressibility of 1^00=230 MeV. 
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FIG. 12: (Color online) Same as in Fig. IIOI but for the cen- 
troids mi /mo. 
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FIG. 13: (Color online) Same as in Fig. 1121 but for the zero- 
range pairing force. 



Bogolyubov solutions and Quasiparticlc-Random-Phase- 
Approximation monopole strength functions in semi- 
magic nuclei. Results were systematically compared with 
those pertaining to the standard zero-range pairing inter- 
action. 

From the monopole strength functions, we extracted 
the finite-nucleus incompressibilities and compared them 
to experimental data. It turned out that neither zero- 
range nor separable pairing effects were able to describe 

TABLE I: Parameters (in MeV) of the liquid-drop formula 
P2p with standard errors, obtained by a fit to the values of Ka 
calculated in M semi-magic nuclei across the mass chart. The 
parameter x was determined as the square root of the sum of 
fit residuals squared divided by the number of fit degrees of 
freedom (M — 5 in our case). 
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the low values of incompressibilities measured in tin, rel- 
ative to the high value measured in ^°^Pb. By changing 
the infinite-matter incompressibility, one can certainly 
describe either the tin or lead values; however, the high 
difference thereof remains unexplained. 

We have also performed adjustments of the liquid-drop 
formula to microscopically calculated incompressibilities, 
and we found that (i) such a formula is able to describe 
microscopic results very well, and (ii) the volume liquid- 
drop term is significantly higher than the infinite-matter 
incompressibility determined for a given functional. 



VI. CONCLUSIONS 

In this work we have presented the first applica- 
tion of the separable, finite-range pairing interaction 
of the Gaussian form together with the non-relativistic 
functional of the Skyrme type. This interaction was 
used to determine both the ground-state Hartrce-Fock- 
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